Spin-Boson Theory for Magnetotransport in Organic Semiconducting Materials 
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We present a spin-boson theory for magnetotransport in organic semiconducting materials, on 
the basis of a coupling between charge carriers' spin and a local bosonic environment, which is 
shown to be an irreducible ingredient in understanding of the anomalous organic magnetoresistance 
(OMR). Among those composing this environment triplet-excitons play a basic role. The incoherent 
hopping rate between molecules is calculated to give out the fundamental behavior of OMR. The 
underlying mechanism is revealed from the calculation of entanglement, represented by the von 
Neumann entropy, between the carrier's spin and bosons. We also obtain the dependence of OMR 
on the bias voltage, the spin-boson coupling, and the boson frequency. The results obtained from 
the theory are in good agreement with experiments. 

PACS numbers: 73.43. Qt, 72.20.Ee, 71.35.-y 



The discovery of room-temperature, low-field mag- 
netoresistance (MR) in organic light-emtting devices 
(OLEDs) was one of the major achievements of spin- 
tronics in the last decade. 0, IH Compared to its inor- 
ganic counterparts, a sizable organic MR (OMR) is rela- 
tively easy to be obtained, showing extensive potential in 
magnetically controlled applications. Q Yet, after years 
of intense research, the origin of this magnetic field effect 
out of these nonmagnetic materials remains controversial. 
Traditional theories of MR in inorganic materials, such 
as that based on Lorentz force or spin-orbit coupling, are 
recognized to be invalid in organic materials. [4j The hy- 
perfine interaction (HFI) caused by hydrogen nuclei in 
organic molecules has been thought to be significant, [5| 
which is verified very recently by an experiment on iso- 
tope effect in spin response of 7r-conjugated polymers. 
Meanwhile, a similar experiment in the Alq3 showed the 
response is isotope independent, 0] which seems interac- 
tions other than HFI may dominate the spin process. 
A bipolaron model[8j using Monte Carlo (MC) simula- 
tion gave a mechanism for both of positive and negative 
OMR at high density of charge carriers. The MC method 
has also been successfully extended to study spin diffu- 
sion in spin valve effect based on incoherent hopping of 
a charge carrier and coherent precession of its spin.Q 
Moreover, there are experiments indicating that the be- 
havior of OMR at low field might be different with that 
at relatively high field. 0, [l(| In a word, a comprehen- 
sive understanding of the experimental phenomena will 
be enlightening. 

A proper theory that describes the intrinsic physics 
of OMR should incorporate the following considerations: 
(i) A prominent characteristic of charge transport in or- 
ganic semiconducting materials is the incoherent hopping 
between molecules; [ll| (ii) The spin decoherence time 
is about 0.5/iS, 12J much longer than the time (~ 100ns 
in Alqa) for carriers transporting in devices. |1 3| Fur- 
ther, the intermolecular hopping was observed to be 
spin-conserving , which implies the motion of carriers' 



spin is coherent; [H (iii) The OMR is robust under room 
temperature. 15[ Even though the Zeeman splitting en- 
ergy of about lOOmT is much smaller than the thermal 
energy, the effect of magnetic field will not be smeared 
out by the thermal fluctuation, for the motion of spins 
being coherent in the whole process. For the same rea- 
son, the effect is little influenced by all other things that 
are not interacting with spins; (iv) There are contradic- 
tions among various experiments on the role of excitons. 
On one hand, these processes associated with excitons, 
such as intersystem crossing between singlet and tri plet 
excitons 0, Ho |. and triplet-triplet annihilation 17 , 18 1 
were regarded as the main response to magnetic field, 
which implies the role is essential. This might be com- 
plemented by a recent experiment of transient electro- 
luminescence, which revealed that OMR could be ob- 
served only when the carriers and excitons coexist in the 
system. [19} On the other hand, there are also evidences 
that the excitons are irrelevant to the OMR, for example, 
the measurement of the intersystem crossing rate shows 
it is magnetic field independent. [3, As a result, it 
can be realized from above that all mentioned mecha- 
nisms overlook the scattering of carriers from excitons, 
which should be a primary role of excitons in the mag- 
netotransport. 

With those considerations, we propose a spin-boson 
theory for magnetotransport in this Letter. The Hamil- 
tonian we suggest for the theory is expressed as 



(1) 



where the intramolecular part of the n-th molecule 

a 

+ gf i B B ■ S„ (2) 

describes the spin interaction of a charge carrier with a 
local bosonic environment in an external magnetic field 
B. Here, S n (S*) is the (z-direction) spin operator of a 
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carrier, b\ a (b n . a ) creates (annihilates) a boson that is 
the a-th mode of the environment, uj n ^ a the correspond- 
ing frequency of bosons, 7„ iQ the spin-boson coupling, 
g the Lande factor that is set to be the commonly ac- 
cepted value 2.0, \ib the Bohr magneton, and z direction 
is chosen along that of the specific bosonic mode at each 
molecule. As an important characteristic of organic ma- 
terials, when a carrier hops into a molecule, it is immers- 
ing in a complex and disordered surrounding medium. 
This medium could be treated as an environment com- 



posed of a number of bosons. 21] Many factors contribute 
to the medium, such as local molecular vibration modes, 
hydrogen nuclear spin, and excitons, in a way that de- 
pends on materials. [6j Especially, triplet excitons should 
play a basic role for their long lifetime and weak diffusive 
ability compared with that of charge carriers. [22J The en- 
vironment could be regarded to be an effective magnetic 
field defined as 
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(3) 



to the carriers' spin plus quantum fluctuation that will 
be shown to be crucial to give the magnetic field effect. 
The intermolecular part of ([T]) 



H' 



,<y + h.C. 



(4) 



describes the carrier's hopping between molecules. 
c li(n') a ( c n(n'),a) creates (annihilates) a carrier with spin 
a at n(n')-th molecule, and J nn i the overlap integral of 
wavefunctions between molecules. This term includes all 
the magnetism independent factors in conventional treat- 
ment of organic charge transport, such as disorders of 
molecular energies and intermolecular distances, which 
ensures that we could only consider a spin-related Hamil- 
tonian for H n . 

Considering the characteristic of organic materials as 
discussed above, we could write the density matrix of 
the system as a direct sum of all local density matrices 
of molecules, that is, 



P = Pl 93 Pi 93 • • • Pn 



(5) 



and divide the whole process into two steps, intramolec- 
ular evolution and intermolecular hopping, to calculate 
the incoherent hopping rate of a single charge carrier. 
The dynamical evolution follows the equation of motion, 



^-Q^Pn = [H„,p n ]. 



(0) 



H n could be numerically diagonalized by constructing 
the following basis, indexed by a set of integer numbers 

{i*}M 

\{l a },8) =]l^ 7 f-(bi + su a ) l -e- su ^\{0}, S ), (7) 



where u a = ja/huJa is the displacement of bosons, 
s = +1 for spin up and —1 for spin down, and for sim- 
plicity we have dropped the molecular index n. Here, u a , 
depending on the spin-boson coupling, is approximately 
equal to the Huang- Rhys factor S (= X/Tiuj with A is the 
reorganization energy), which is commonly of the order 
~ 1 in organic materials. [2~ij Throughout this work, we 
set the cutoff number of bosons to be 80 to ensure the 
calculation is convergent. 

The second step is the incoherent hopping between 
molecules, whose rate could be derived in terms of the 
Fermi golden rule. (25j 



v*^\(f\H'\i) 



(8) 



where r if = Ah 2 sm 2 [(E f - E l )t d /2h]/(E f - Ei) 2 t d with 
td being the decoherence time within which the hopping 
is coherent, and pi and pf (|/)) are the initial and 
final density matrix (state) with the energy Ei and Ef, 
respectively, expressed in Eq. (|5|). td is a quantity de- 
termined approximately by the molecular structure re- 
laxation, which gives it's at the order of O.lps. 11( It's 
noted that this treatment is qualitatively similar to the 
Franck-Condon principle, 21| which has been widely used 
in organic electronics as well as the spin-dependent ex- 
citon formation. [2^ To determine the initial and final 
state in ([8]), we consider the following process. At a time 
when the process begins, a carrier with some spin s hops 
onto the n-th molecule, and the density matrix could be 
expressed as a direct product 



Pn(t = 0) = g n (s) ® g„ 



(9) 



where g denotes the (2x2) spin density matrix of a car- 
rier while g for bosons. Right at this moment, the in- 
teraction between the carrier and local bosonic modes is 
switched on, and then they become entangled. The sys- 
tem evolutes following Eq. ([6]) till the carrier hops to 
the unoccupied n'-th molecule at time t w , called waiting 
time, whose average value could be estimated by the mo- 
bility of carriers and bias voltage. For organic materials, 
t w 3> td- After the hop, the carrier's spin and bosons 
on the n-th molecule disentangle. The spin state keeps 
unchanged because of spin-conserving, while the bosons 
will reorganize. An adiabatic elimination procedure tech- 
nique is used for the spin state, 27] which gives the spin 
density matrix at the n'-th molecule as 



Qn'(s) = Tr b {Pn(t w )}, 



(10) 



where the trace with subscript b is that over all bosonic 
degrees of freedom. For the reorganization of bosons, it 
depends on the relative ratio of the relaxation time of the 
bosonic environment (t r ), which is not explicitly known, 
and that of molecular structure (~ t d )- In case t r 3> td, 
the relaxation is negligible, so the bosonic density matrix 
becomes 



Qn 



Tr s {/) n (t w )}, 



(11) 
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FIG. 1: Magnetic-field dependence of (a) the hopping rate v 
(in unit of rj 2 /h 2 ) between two molecules and (b) the von 
Neumann entropy for various waiting times t w . 7 = l.O^ieV 
and huj = 2.0fieV. 
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FIG. 2: (a) Magnetic field dependence of OMR under various 
bias voltages, (b) The dependence of OMR on the voltage 
under various magnetic fields. Inset shows the OMR is insen- 
sitive to 77. 



where the trace with subscript s is that over spin degree 
of freedom. But for t r <C td, the relaxation is complete, 
and g n should return to its equilibrium state, in which 
all u a = 0. For an arbitrary ratio of t r /td, we could 
introduce phenomenological parameters rj a depending on 
cu a between and 1. Thus, g n in Eq. ([7]) should be that 
with all u a replaced by m*(e r) a u a ). It will be seen later 
that the result is insensitive to the choice of r\ a , i.e., the 
relaxation of bosonic environments. 

Now, the hopping rate could be derived directly from 
Eq. (O as 

v^n- = Tnn '/ 2 nn ' r Tr{g nQn ,(s)p n (t w )}, (12) 
n 

since it concerns only two molecules. Considering the 
random directions of a carrier's initial spin and its cou- 
pling with bosons, the hopping rate we calculate below 
will be averaged over all directions of the carrier's spin 
and that of applied magnetic field since the spin-boson 
coupling has been fixed along the z-direction. 

Within the framework of perturbation theory, the in- 
teraction between different bosonic modes is not signifi- 
cant, so it is meaningful to investigate the system with an 
environment of a single bosonic mode. Hereafter we omit 
those molecular and boson's indexes if it would not create 
confusion. At first, the calculated hopping rate is shown 
in Fig. [Ua), from which we can see its strong depen- 
dence on external magnetic field as expected. We have 
set 7 = 1.0/xeV and hoo = 2.0/ieV, which correspond to 
an effective magnetic field Bq(= 2j 2 /hwgfiB) — 8.6mT 
evaluated from Eq. (3). This field strength is similar 
with that used in HFI model for OMR[5J and extracted 
from muon experiment [28]. Accordingly, a saturation is 
reached for the magnetic field being 50 - lOOmT depend- 
ing on the waiting time t w . It is clearly seen that the 
magnetic field dependence becomes stronger with a larger 
t w . The slow waviness in those curves is due to the Rabi 
oscillation between spin states. 

The observed magnetic field dependence originates 
from the incoherent hopping, which disentangle the 
spin and bosonic environment. To show it we calcu- 
lated the von Neumann entropy [2^], which is defined as 



Tr(f51og 2 g). The result in Fig. QJb) implies that the en- 
tanglement weakens the ability of a carrier hopping out 
of the environment. The external magnetic field changes 
the entanglement between the carrier and the bosonic en- 
vironment, and then the magnetic field effect arises. This 
is the basic mechanism of OMR from our theory. 

With the hopping rate between molecules in Eq. (12), 
we have the OMR as 

AI 1(B) 1(0) (u(B)) MO)) m , 
/ - 1(0) (K0)) ' 1 ' 

where 1(B) is the current depending on the magnetic 
field B, and (•) is the average over samplings with a 
distribution of t w here. [2^] The average of t w should be 
proportional to the inverse of bias voltage applied to a 
given device, which is expressed as (t w ) = Lj^l/pV^, 
with Lj) being the thickness of device, I the intermolec- 
ular distance, fi the mobility of carriers, Vj, the bias volt- 
age, and £ a factor introduced to account for the car- 
rier's random hop. For example, we consider a typical 
OLED structure ITO (indium tin oxide)/NPB(7V, AT'-di- 
1-naphthyl-iV, iV'-diphenylbenzidine) (50 nm) / Alq3 (tri- 
(8-hydroxyquinolinc)-aluminum) (50nm)/ Al (aluminum). 
The effective mobility of electrons in Alq3 is taken as n s» 
10 _4 cm 2 /Vs, the average distance between molecules I « 
0.5nm, and £ « 5, so we have, e.g., (t w ) s» 0.25ns for Vb = 
4.0V. Fig.[5[a) shows the OMR under different bias volt- 
ages. The basic line shape obtained here is in good agree- 
ment with a consensus based on experiments, fll. IBl fl9L l3pj 
Furthermore, we show the dependence of OMR on bias 
voltage under different magnetic field in Fig.[2][b). A de- 
cay following the voltage increase is found, which matches 
experimental observations, [l^] For a small bias voltage, 
e.g., \% < 2.0V, the current is extremely small and the 
scattering between excitons becomes important, which 
is not included in this work, that might be the origin 
of positive OMR observed in experiments. |3G] The rj de- 
pendence of OMR is given in the inset, which verifies the 
choice of rj has little influence. 

Now we come to see the influence of a bosonic mode 
on the OMR. Fig. [3][a) shows the dependence of satu- 
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FIG. 3: OMR dependence on (a) the effective magnetic field 
Bo and (b) the boson energy hw with t w = 0.25ns. 



OMR and its dependence on bias voltage in consistent 
with a consensus based on experiments. Triplet excitons 
are believed to contribute mostly to this environment, 
and relevant experiments are understood in accordance 
with the theory. 
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sions. The work was supported by the NSF of China, the 
National Basic Research Program of China (Grant No. 
2009CB929204), and the EC Project OFSPIN (Grant No. 
NMP3-CT-2006-033370) . 



rated AI /I on the spin-boson coupling 7, which is equiv- 
alent to the effective magnetic field Bq, at t w — 0.25ns. 
Obviously, OMR is found to be nearly proportional to 
Bq in the range we calculated. Fig. EJb) shows the 
frequency dependence of OMR. As expected, OMR ap- 
proaches to vanishing when the frequency becomes large 
enough, where the bosonic mode is equivalent to an ef- 
fective mag netic field, that is, the (classical) hypcrfinc 
field. [1, Il0| From the perspective of entanglement, it's 
easy to understand the vanishing OMR since a classi- 
cal field does not entangle with the carriers' spin. This 
also shows the quantum nature of local bosonic environ- 
ments is irreducible for OMR.[3l| It appears that the 
OMR peaks all reside at Tiuj 5.0/xeV, which is caused 
by setting the waiting time at a fixed value t w = 0.25ns. 
Such phenomenon implies that bosonic modes of com- 
patible energy with the waiting time contribute more to 
OMR. 

Until now, we have discussed the environment of a sin- 
gle bosonic mode. In reality, a number of modes should 
contribute to OMR, and then the effective magnetic field 
B we used in the one-mode case could be regarded as 
that defined in Eq. (3) within the perturbation theory. 
As has discussed, triplet-excitons play a basic role in the 
environment, so Bq should be strongly dependent on the 
density of excitons that could be adjusted experimentally 
in many ways. Moreover, another important physical 
quantity to OMR is the boson energy %u>, which should 
depend on the constitution of the environment, related 
to characters of specific material, such as the mass of 
atoms and/or radicals in the molecule, electronic excita- 
tion energies. Due to the nonmonotonicity of OMR on Tiuj 
as shown in Fig. |3ljb), the discrepancy of isotope effect 
in different materials [fl 0] might be understood within 
the theory. It is also shown[l7( experimentally that the 
OMR is strongly associated with the excitation energy 
by changing dopants. 

In summary, we have proposed a theory based on the 
coupling between carriers' spin and a local bosonic en- 
vironment. The entanglement is calculated to reveal 
the underlying mechanism of the anomalous OMR. The 
quantum nature of the environment is shown to be irre- 
ducible. The theory establishes the basic line shape of 
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